SEASONAL CHANGES IN TITAN'S SURFACE TEMPERATURES 
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ABSTRACT 

Seasonal changes in Titan's surface brightness temperatures have been observed by 
Cassini in the thermal infrared. The Composite Infrared Spectrometer (CIRS) 
measured surface radiances at 19 Jim in two time periods: one in late northern 
winter (L s = 335°] and another centered on northern spring equinox ( L s = 0°]. In 
both periods we constructed pole-to-pole maps of zonally averaged brightness 
temperatures corrected for effects of the atmosphere. Between late northern winter 
and northern spring equinox a shift occurred in the temperature distribution, 
characterized by a warming of ~0.5 K in the north and a cooling by about the same 
amount in the south. At equinox the polar surface temperatures were both near 91 K 
and the equator was 93.4 K. We measured a seasonal lag of A 9° in the 
meridional surface temperature distribution, consistent with the post-equinox 
results of Voyager 1 as well as with predictions from general circulation modeling. A 
slightly elevated temperature is observed at 65° S in the relatively cloud-free zone 
between the mid-latitude and southern cloud regions. 


Key words: — infrared: solar system — planets and satellites: individual [Titan] — 
radiative transfer — radiation mechanisms: thermal 


Goddard Space Flight Center, Greenbelt, MD 20771; donald.e.jennings@nasa.gov 

2 University of Maryland, College Park, MD 20742 

3 The Catholic University of America, Washington, DC 20064, USA 

4 LESIA, Observatoire de Paris-Meudon, 92195 Meudon Cedex, France 

5 Universitat zu Koln, Albertus-Magnus-Platz, 50923 Koln, Germany 


1. INTRODUCTION 

A primary goal of the Cassini mission is to study changes on Titan as it 

progresses from northern winter to summer. Major variations in Titan's weather 

can be expected during the extended mission, driven by seasonal shifts in surface 

heating. The small amount of sunlight that makes it through the atmospheric haze 
warms the surface, which in turn heats the lower troposphere and gives rise to 




atmospheric convection and circulation. Surface warmth lifts volatile materials, 
which condense to produce clouds and rain (Griffith etal 2000, Brown etal. 2002). 
Meridional distributions of precipitation evolve globally as the insolation cycles 
through the Saturnian year. Modified appearances of surface features observed by 
Cassini have been associated with seasonal storms (Turtle etal. 2011, Turtle etal 
2009), and the appearance of southern clouds has diminished in both Cassini and 
ground observations as Titan has progressed into southern autumn (Schaller etal 
2006, Rodriquez et al 2009). In addition, warming of the surface governs wind 
patterns including dune-forming wind reversals and a seasonal shift of the ITCZ 
(intertropical convergence zone) that occur at the equinoxes (Mitchell et al. 2006, 
2009, Tokano 2010). Seasonal changes in the atmosphere and on the surface are 
likely to accelerate as Titan moves toward northern summer. 

Voyager I first measured the temperature of the surface of Titan through a 
spectral window of low atmospheric opacity at 19 pm (Samuelson etal 1981, Flasar 
etal 1981). From the single Voyager I encounter it was not possible to derive 
information about seasonal change, but the measurements showed a symmetric 
distribution around the equator as would be expected for the time of the flyby just 
after northern spring equinox. The Composite Infrared Spectrometer (CIRS; Flasar 
etal 2004) aboard Cassini covers the 19-pm spectral window and has been 
recording spectra of Titan since 2004. We previously reported CIRS results from the 
2004-08 portion of the mission (Jennings etal 2009) that gave a temperature of 
93.7 K at the equator, in agreement with Huygens (Fulchignoni et al 2005). The 
temperatures were 2 K cooler at the south pole and 3 K cooler at the north pole, 
appropriate to late northern winter. Titan has since advanced through northern 
spring equinox and into early northern spring, and the meridional distribution of 
surface temperatures is expected to have shifted northward significantly during that 
time. We present here evidence that the surface temperature is responding to the 
seasonal insolation shift and that the changes follow predictions from general 
circulation models (GCM; Tokano 2005). Recently, Cottini etal (2011) have 
reported spatial and temporal variations of surface temperature from CIRS 
measurements, including evidence of diurnal variations and seasonal shifts. Here we 
extend their investigation of seasonal changes. 

2. OBSERVATIONS 

Surface brightness temperatures were derived from far-infrared observations by 
CIRS following the method described by Jennings et al. (2009). Spacecraft distances 
were constrained to be within 140,000 km of Titan to keep the 3.5 mrad field of 
view smaller than about 10% of the disk. The average emission angle viewed was 
restricted to <70° to limit the effect of the atmosphere. The surface radiances were 
measured in a channel centered on 530 cm ' 1 with a bandwidth set by the 15 cm ’ 1 
apodized resolution of the instrument. This spectral position is well within the 19- 
micron atmospheric window between the CH4-N2 and H 2 -N 2 collisionally induced 
opacities (Courtin etal 1995; Samuelson etal 1997) and separated from HC 3 N and 



C 2 HD molecular emission below 520 cm' 1 (Kunde etal, 1981; Coustenis etal, 

2008). The choice of 530 cm' 1 also permits direct comparison with previous work 
(Flasar etal 1981, Jennings etal 2009). On each Titan flyby only a limited latitude 
interval is observed and more than a year is typically required to cover the full pole- 
to-pole range. We have identified two distinct periods that each contain sufficient 
data to sample all latitudes and which are separated sufficiently in time to show a 
seasonal shift. These two periods are September 2006 through May 2008, hereafter 
referred to as late northern winter (LNW), and November 2008 through May 2010, 
hereafter called northern spring equinox (NSE). The periods are centered at solar 
longitudes L s - 335° (0.8 Saturn month before NSE) and L s = 0° (at NSE). 

3. DATA ANALYSIS 

Radiance spectra in each time period were averaged zonally in 10° latitude bins 
before being converted to brightness temperature. We used larger latitude bins than 
the 5° bins in our previous study (Jennings etal 2009) to improve the sensitivity of 
the results to seasonal change. Spectra in each latitude bin were averaged in two 
ranges of emission angle, 0-50° and 50-70°, and these were corrected separately for 
the effects the atmosphere. To correct for the atmosphere we calculated the opacity 
and emission from the surface to space using the method described by Jennings et 
al (2009), with several modifications to update our model atmosphere. As in the 
earlier work, we adopted the Huygens Atmospheric Structure Instrument (HASI) 
temperature profile for the equatorial region (Fulchignoni etal 2005). This profile 
covers 0-147 km altitude. We then modified the temperature profile for higher 
latitudes to match the altitude behavior reported from Cassini Radio Science 
Subsystem (RSS) occultations (Flasar and Schinder 2010) and from CIRS molecular 
line retrievals (Achterberg etal. 2008, Achterberg etal 2011, Anderson etal 2010, 
Coustenis et al 2010). We now include in our temperature model: 1) a large 
stratospheric variation in the north with a local maximum near 90 km and a local 
minimum near 120 km; 2) a middle troposphere (10-30 km) independent of 
latitude; and 3) a layer below 10 km where the lower atmospheric temperature 
transitions to the surface temperature. We adopted a haze opacity that increases 
downward with a scale height of 65 km above 80 km and remains constant from 80 
km to the surface (Cottini etal 2011, de Kok etal 2007, Tomasko etal 2005). In 
addition, we applied a latitude-dependent haze factor increasing by 50% from the 
south to the north, with a maximum at 40N (Cottini et al. 2011, Rannou et al 2010, 
Samuelson et al 1997). The haze opacity was adjusted by an overall factor to 
produce the dependence of radiance on emission angle that we measure at the 
equator and at 73° N (Jennings et al 2009), In our calculation of CH 4 -N 2 opacity the 
methane mole fraction was updated to the dependence reported by Niemann etal 
(2010): 1.48% above about 40 km rising to 5.65% below 7 km. The temperature 
dependent absorption coefficients for CH 4 -N 2 from Borysow & Tang (1993) were 
increased by 50% as recommended by de Kok, Irwin &Teanby (2010). For the H 2 - 
N 2 opacity we used a mole fraction for molecular hydrogen (H 2 ) of 0.001 (Courtin et 
al. 1995, Jennings etal 2009, Niemann etal. 2010) and assumed that it was constant 
over altitude and latitude. Temperature-dependent absorption coefficients for H 2 -N 2 



were taken from Courtin [1988] and Dore etal. (1986). We assumed the emissivity 
of the surface at 19 microns to be unity, which is justified because: 1) our measured 
surface brightness temperature near the equator is the same as that measured at the 
surface by Huygens (93.7 K); and 2) our results are in agreement with near-surface 
radio occultation measurements (Flasar and Schinder 2010). 

4. RESULTS AND DISCUSSION 

Figure 1 shows our derived surface brightness temperatures for the two periods. 
The temperatures for LNW are similar to those we previously reported for the first 
portion of the Cassini mission (Jennings etal 2009). At the equator the temperature 
was -93.6 K, while at the south and north poles it was -91.5 and -90.5 K, 
respectively. However, between LNW and NSE it is evident that the temperatures 
increased in the north by roughly 0.5 K and decreased in the south by about the 
same amount. Equatorial temperatures may have decreased slightly, by about 0.2 K, 
possibly because Titan moved away from perihelion during this period. Although 
the maximum surface temperatures for both periods were near the equator, both 
profiles had generally higher temperatures in the south. The centroid of the 
temperature distribution, determined from cosine fits (as in Jennings etal. 2009), 
was located at 15° S latitude for LNW and 4° S for NSE. This gives a seasonal lag 
behind insolation of ALs- 9°. Our surface brightness temperatures are similar, in 
both seasonal lag and magnitude, to those of Voyager 1 IRIS (Flasar et al 1981, 
Courtin & Kim 2002). Voyager 1 observed a roughly symmetric distribution about 
the equator at Ls= 10° with temperatures ranging from about 91 to 93 K for ±60° 
latitude. We conclude that the surface temperatures are repeating after one Saturn 
year, both in amplitude and seasonal lag. 

Figure 2 shows our surface brightness temperatures compared with two surface 
composition scenarios of Tokano (2005): porous ice regolith and rock-ice mixture. 
Both LNW and NSE are compared with the GCM predictions for the same solar 
longitude. We find that the distribution of the observed surface temperatures has 
tracked the predictions quite closely. The distinction between the two scenarios lies 
in the thermal inertia, which is eight times higher for rock-ice mixture than for 
porous ice regolith (with an assumed thermal inertia of -330 J m~ 2 s" 0 ' 5 K~') due to 
the high heat conductivity of ice and higher compaction. In both time periods the 
measured equator-to-pole variations match the porous ice regolith scenario more 
closely than the rock-ice mixture. However, while the observations are closer to 
porous ice regolith in the south, in the north they tend to lie between the two 
scenarios. Although this might imply a difference in composition between the north 
and south, we are not sure that the accuracies of the data or the model warrant such 
a conclusion. In terms of seasonal lag, a cosine fit to the GCM predictions for a 
porous ice regolith shows a north-south symmetry at Ls~ 13°, in agreement with 
our derived A Ls~ 9° within uncertainties. In contrast, the predictions for a rock-ice 
mixture are symmetric at Ls~ 30°, inconsistent with both our results and those of 
Voyager 1 IRIS. 



The temperatures at the poles in both the LNW and NSE measurements are near 
the triple points for methane (90.6 K) and ethane (89.9 K), suggesting that gaseous 
and condensed phases might exist together on the surface at high latitudes. A ~0.5 K 
change (our observed increase in the north and decrease in the south) corresponds 
to a 7% change in saturation vapor pressure for methane and 14% change for 
ethane. Evaporation and condensation rates are therefore likely to have already 
been affected by the seasonal advancement. By the end of the Cassini mission, at 
summer solstice, the poles may have changed by as much as 3 K from the LNW 
values (Tokano 2005). This represents a 50% change in vapor pressure for methane 
and a factor of 2 for ethane. Thus, we can expect major changes in condensed phases 
at high latitudes during the remainder of the mission. 

Some features in our observed surface temperatures may be related to cloud 
distribution. In the south the latitude bin at 60-70° S is elevated in temperature 
relative to its neighbors during both LNW and NSE, by about 0.5 and 0.7 K, 
respectively. We note that this latitude coincides with a relatively cloud-free gap 
between two zones where clouds commonly appear, one zone centered at ~40° S 
and another region extending poleward from -70° S (Rodriquez etal. 2009, Griffith 
et al. 2005, Roe et al. 2005). An elevated surface temperature at 65° S might be the 
result of higher insolation in the cloud-free zone. As northern summer approaches 
the atmospheric circulation patterns are expected to change and the temperature 
peak at 65° S may migrate in latitude or disappear. 

In the north the largest temperature changes are at higher latitudes, >40° N. In 
fact, we see no changes greater than the errors in the two bins centered at 25 and 
35° N, possibly because the increase in solar declination is compensated by the 
increasing solar distance. The temperature variation is greatest at 50-70° N, but is 
small again at 70-90° N. At these highest latitudes the changes are less than 
predicted for land composed only of porous ice regolith. This may result in part 
from the large fraction of lake surface at northernmost latitudes. During the winter 
the polar lakes are expected to have a temperature similar to that of the land. 
However, in the spring the warming of the lakes is slower than the land due to their 
larger thermal inertia (Tokano 2009). The predicted lake temperature is 90-91 K 
throughout winter and early spring and begins to increase only in late spring. Our 
results so far are consistent with this prediction, implying that the large lake 
coverage effectively increases the thermal inertia near the pole. By the end of the 
Cassini mission we may begin to see evidence of rising temperatures in the polar 
lakes. 

We acknowledge support from NASA’s Cassini mission and Cassini Data Analysis 
Program. TT was supported by DFG. 


REFERENCES 



Achterberg, R. K., Conrath, B. J., Gierasch, P. J., Flasar, F. M., & Nixon, C. A. 2008, 
Icarus, 194, 263 

Achterberg, R. K., Gierasch, P. j., Conrath, B. ]., Flasar, F. M., & Nixon, C. A. 2011, 
Icarus, 211, 686 

Anderson, C.M., Samuelson, R.E., Bjoraker G.L., Achterberg, R.K. 2010, Icarus, 207, 
914 

Brown, M. E., Bouchez, A. H„ and Griffith, C. A. 2002, Nature, 420, 795. 

Borysow, A., & Tang, C. 1993, Icarus, 105, 175 
Courtin, R. 1988, Icarus, 75, 245 

Courtin, R., Gautier, D., & McKay, C P. 1995, Icarus, 114, 144 
Courtin, R., & Kim, S. J. 2002, Planet and Space Sci., 50, 309 
Coustenis, A. etal, 2008, Icarus, 197, 539 
Coustenis, A. etal 2010, Icarus, 207, 461 

Cottini, V., Nixon, C. A., Jennings, D. E., de Kok, R., Teanby, N. A., Irwin, P. G., & Flasar, 
F. M. 2011, Planetary and Space Science, in press. 

Dore, P., Borysow, A., & Frommhold, L. 1986, J. Chem. Phys., 84, 5211 
de Kok, R., Irwin, P.G.J., Teanby, N.A, Nixon, C.A., Jennings, D.E., Fletcher, L, Howett, 
C., Calcutt, S.B., Bowles, N.E., Flasar, F.M., & Taylor, F.W. 2007, Icarus. 191, 
223. 

de Kok, R., Irwin, P.G.J., & Teanby, N.A 2010, Icarus, 209, 854 
Fulchignoni, M., etal 2005, Nature 438, 785 

Flasar, F. M., Samuelson, R. E., & Conrath, B. J. 1981 Nature, 292, 693 
Flasar, F. M„ etal 2004, Space Sci. Rev., 115, 169 

Flasar, F. M., and Schinder, P. J. 2010, “The thermal structure of Titan's atmosphere 
from Cassini radio occupations” in Through Time; A Workshop On Titan's 
Past, Present and Future, at NASA Goddard Space Flight Center, April 6th - 
8th, 2010. Edited by V. Cottini, C. Nixon, and R. Lorenz., p.51 
Griffith, C. A., Hall, J. L., & Geballe, T. R. 2000, Science, 290, 509 
Griffith, C. A. etal 2005, Science, 310, 474 
Jennings et al. 2009, Astrophys J. Lett., 691, L103 
Kunde, V. G., Aikin, A. C., Hanel, R. A., Jennings, D. E., Maguire, W. C., & 

Samuelson, R. E. 1 981, Nature, 292, 686 

Mitchell J. L., Pierrehumbert, R. T, Frierson, D. M.W., Caballero, R. 2006, PNAS, 103, 
18421 

Mitchell J. L., Pierrehumbert, R. T, Frierson, D. M.W., Caballero, R. 2009, Icarus, 203, 
250 

Niemann, H. B., Atreya, S. K., Demick, J. E., Gautier, D., Haberman, J. A., 

Harpold, D. N. ; Kasprzak, W. T., Lunine, J. I. , Owen, T. C, & Raulin, F. 2010, J. 
Geophys. Res., 115, E12006 

Rannou, P., Cours, T., Le Mouelic, S., Rodriguez, S., Sotin, C., Drossart, P., Brown, R. 

2010, Icarus, 208, 850 
Rodriquez etal. 2009, Nature, 459, 687 

Roe, H. G., Brown, M. E., Schaller, E. L., Bouchez, A. H„ Trujillo, C. A. 2005, 

Science, 310, 477 

Samuelson, R. E., Hanel, R. A., Kunde, V. G,, & Maguire, W. C. 1981, Nature, 292, 688 
Samuelson, R. E., Nath, N. R., & A. Borysow 1997, Planet. Space Sci., 45, 959 



Schaller, E. L., Brown, M. E. ; Roe, H. G., Bouchez, A. H. & Trujillo, C. A. 2006, 
Icarus, 184, 517 

Tokano, T. 2005, Icarus, 173, 222 
Tokano, T. 2009, Icarus, 204, 619 
Tokano, T. 2010, Aeolian Res. 2, 113 

Turtle, E. P., Perry, J. E., McEwen, A. S., DelGenio, A. D., Barbara, j. M., West, R. A., 
Dawson, D. D., and Porco, C. C. 2009, Geophys, Res. Lett. 36, L02204 
Turtle, E. P., Perry, J. E., Hayes, A. G., Lorenz, R. D., Barnes, J. W., McEwen, A. S., 
West, R. A., Del Genio, A. D,, Barbara, {. M., Lunine, j. I., Schaller, E. L., Ray, T. 
L., Lopes, R. M. C., Stofan, E. R. 2011, Science 331. 1414 


Figure Captions 

Figure 1. Titan surface brightness temperatures measured in late northern winter 
(LNW, Sep 2006 to May 2008) and around northern spring equinox (NSE, Nov 2008 
to May 2010). Temperatures were zonally averaged in 10° latitude bins and the 
error in each bin is the standard deviation of its average. The latitude of each data 
point is the average latitude for data in that bin. Southern latitudes are negative. 

Figure 2. Measured surface brightness temperature distributions (blue) for two 
seasonal periods compared with GCM predictions (Tokano 2005) for two types of 
surface material; porous ice regolith (orange) and rock-ice mixture (green). CIRS 
measurements are the same as in Figure 1. a) Late northern winter (LNW, Sep 2006 
to May 2008). b) Northern spring equinox (NSE, Nov 2008 to May 2010). 
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